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The results of tests made with a two dimensional diffusor 
operating with transverse grooves cut into the diverging side walls 
were compared with the results of similar tests made using a smooth 
wall diffusor. The results of this comparison revealed thats 

1 , The efficiency of the ribbed wall configuration exceeded that 
of the smooth wall diffusor for total opening angles of approximately 
28 to 33 degrees. 

2 . The ribbed wall configuration had an efficiency approximately 
120 per cent of that of the smooth wall diffusor when both were oper- 
ated at a total opening angle of 32 degrees. 

3 » The efficiencies of the ribbed wall diffusor operating at 
large opening angles varied as a function of the number of grooves in 
the side walls. 

4 . The efficiency increased as the number of grooves increased, 
but an unexplained drop in efficiency occurred at openiTvg angles of 
30 and 32 degrees with an eight groove configuration. 

5. The maximum number of grooves tested was ten. There is no 
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evidence to indicate that this is an optimum configuration. 
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CHAPTER I 



INTRODUCTION 

Much work has been done and many papers written concerning studies 
of diffusor performance* Most of these studies are limited to the 
study of diffusors operating with a total opening angle of less than 
1^ degrees. Cocanower, Kline, and Johnston, Ref. S have proposed a 
general theory for predicting diffusor performance for a variety of 
two and three dimensional diffusors. Their theory has been shown to 
be in good agreement with experimental tests. 

The studies of diffusor performance for diffusors operating at 
large opening angles have not progressed as rapidly as the small angle 
diffusor work because of the complications associated with the severe 
turbulence resulting when flow separation occurs at the diffusor walls. 
This separation of flow is responsible for the large losses and poor 
performance normally associated with wide angle diffusors. 

Two studies which deal with the problem of wide angle diffusors 
are the papers of Waitman, Reneau, and Kline, Ref. 2 , and of Moore, 
Carl, and Kline, Ref. 3* Each of these studies deals with the condi- 
tions affecting flow separation in wide angle diffusors. 

The flow separation phenomena observed and discussed by Moore, 
Carl, and Kline in connection with water table tests were also observed 
in the air flow tests of this investigation. These phenomena are dis- 
cussed in Section 3*6 of Chapter III. 

Migay, Ref. 4, noted an improvement in the efficiencies of dif- 
fusors operating at wide opening angles when transverse grooves were 
cut into the diverging walls. His explanation for this improvement in 



11 



efficiency was that the "large separation" normally occurring without 
the grooves is replaced, because of their effect, by a series of 
"smaller separations." The resulting losses associated with these 
"smaller separations" are apparently less than the loss associated with 
the "large separation." The efficiency is thereby improved. 

Migay cites the results of many tests which support his claim of 
an improved efficiency when operating with a grooved configuration. 
Unfortunately, many of the results cited are not adequately documented. 
The results of the two dimensional experimentation, for example, are 
presented as the ratio of ribbed wall diffusor efficiency to that of a 
smooth wall diffusor operating under the same conditions. No informa- 
tion is available concerning the specific conditions under which the 
tests were conducted, and the actual efficiencies of each of the dif- 
fusors are not presented in the paper. 

The over-all objective of this study was to determine experimen- 
tally the performance of a two dimensional ribbed wall diffusor. It 
was hoped that the optimum operating range of such a diffusor could be 
found and compared with a smooth walled diffusor operating under simi- 
lar conditions. 

Figure 1 illustrates the type of diffusor used for the experiment. 
The detailed insert in the figure shows what is meant by the term 
"ribbed wall." 

Because of time limitations, the scope of this study was narrowed 
to the study of one particular rib configuration. The geometry of 
this configuration is discussed in Section 2.4 of Chapter II. A limited 
investigation was also made to determine the effect that varying the 
number of grooves had upon the performance of the ribbed wall diffusor. 
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ITEM A 




TWO DIMENSIONAL DIFFUSOR 




FIGURE 1 

TWO DIMENSIONAL DIFFUSOR AND 
DETAIL OF A RIBBED WALL 
CONFIGURATION 
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Migay cites a limited amount of data related to a two dimensional 
ribbed wall diffusor. The data were obtained using a one sided dif- 
fusor with an asymmetric opening angle of 20 degrees. Figure 2 shows 
the type of diffusor used. Table I is a summary of the information 
presented in Ref, 4 related to two dimensional diffusors. 

The initial rib configuration utilized for this study was chosen 
to be approximately equal to that used in Ref, 1, In this study, how- 
ever, a symmetric diffusor was used. 
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ILLUSTRATION OF AN ASYMMETRIC TWO DIMENSIONAL DIFFUSOR OF THE TYPE USED IN REF. 
AND A SYMMETRIC TWO DIMENSIONAL DIFFUSOR OF THE TYPE USED IN THIS STUDY 



TABLE I 



SUMMARY OF TWO DIMENSIONAL DIFFUSOR STUDIES FROM REF. 1 



Effect of Parameter A (B = 5mm) 





1 .0 


0.857 


0.915 


1.25 


1.46 


1.64 


1.68 


1.71 




A (mm) 


0.0 


0.500 
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2.00 


4.00 
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10.0 


15.0 




A/B 


0.0 


0.100 


0.250 


0.400 


0.800 


1 .600 


2.00 


3.00 





Effect of Parameter B (A = 20mm) 





1.61 


1.7 


1.71 




0.5 




B 


1.0 


3.0 


5.0 


unstab. 
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A/B 


20 


6.7 


4.0 




1.25 





Effect of Number of Grooves N (B = 1mm, A = 20mm) 





1.04 


1.13 


1.30 


1.515 


1.63 


1.68 


1.71 


1.71 


1.71 


N 


1 


2 


3 


4 


5 


6 


7 


10 


15 



Effect of Number of Grooves N (B = 3mm, A = 20mm) 





1 .06 


1.14 


1.29 


1 .40 


1.54 


1 .57 


1.63 


1.64 


1.64 


N 
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3 


5 


7 


9 


1 1 


13 


15 


40 



Effect of Number of Grooves N (B = 5mm, A = 20mm) 





0.6 


0.7 


0.78 


1 .5 


1.68 


1.71 


1.71 






N 


1 


2 


3 


4 


5 


6 


15 
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CHAPTER II 



EXPERIMENTAL PROCEDURES 

2*1 Experi menta 1 set -up. The equipment used for this study was 
located in the "New Addi t i on" of Building 230 at the Jet Propulsion 
Laboratories, Naval Postgraduate School, Monterey, California. Figure 
3 is a schematic of the installation. The source of air for the tests 
was the Carrier Air Compressor installed in Building 230, 

P>"^ssure readings . All pressure measurements in the flow 
field were made using a Fused Quartz Precision Pressure cJltge manufac- 
tured by Texas Instruments Inc, of Houston, Texas, Prior to, and twice 
during the time period of the tests, the pressure instrument was cali- 
brated using a water manometer. The calibration remained invariant 
throughout the duration of the test; and the value of 0,8 inches of 
water pressure per increment, or count, on the instrument was used as 
the conversion coefficient for all pressure measurements. The details 
of the conversion of the pressure and temperature readings are dis- 
cussed in Chapter III, Sections 3*1 and 3*2, 

Atmospheric pressures were read from a standard wall barometer. 

Two pressure readings were required at the flow nozzle. The static 
pressure upstream, P^j, and the drop in static pressure across the 
nozzle, AP , were measured and recorded for the purpose of determining 
the flow rate through the system. The computations involved in the 
determination of the flow rate are discussed in Chapter III, Section 

3.3. 

Two pressures were measured at the inlet plane of the diffusor, 
the static pressure and the dynamic head, ^ , The static pressure 
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was assumed to be uniform in the plane of the throat. Two static taps 
were located at the top of the diffusor inlet and two at the bottom. 
All four of these static pressure taps were connected to the pressure 
measuring instrument through a common lead. The total pressure was 
measured at three positions in the central region of the flow, free 
of the boundary layer effect. Figure 4 shows the location of the pres 
sure orifices and the details of the total pressure probe construction 
The probe was manufactured by Mr. M. 0 dell, machinist. Aeronautics 
Department, NPGS. 

Figure 5 is a schematic of the pressure instrument hook-up. By 
opening and closing various pairs of the six switches, it was possible 
to measure all of the required pressures. 

Section 3-4 of Chapter III illustrates how the conditions at the 
throat were calculated using these pressure measurements. 

2.3 Temperature measurements . Two total temperature readings 
were taken in the flow field. The first was taken just ahead of the 
flow nozzle and was used in the flow rate calculation. The second 
total temperature measurement was taken just ahead of the diffusor 
inlet. This temperature was used in the determination of the flow 
conditions at the inlet. 

All flow field temperatures were measured using a thermocouple. 
Section 3.1 of Chapter III discusses the details of the temperature 
reading conversions. The ambient air temperatures were measured with 
ordinary mercury bulb thermometers. 

^^Qsur i ng and recordi ng the data . The rate of flow through 
the system was regulated by bleeding the compressor output to the 
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DIFFUSOR INLET PLANE 




CROSS SECTION TOTAL PRESSURE ORIFICE 
FIGURE 4 

SCHEMATIC OF THE DIFFUSOR INLET PLANE ILLUSTRATING THE LOCATION 
OF THE TOTAL AND STATIC PRESSURE TAPS 
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FIGURE 5 

SCHEMATIC OF THE PRESSURE INSTRUMENTATION INSTALLATION 
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atmosphere. The maximum obtainable flow rate was approximately 2.4 
Ib^sec. At this flow rate, the velocity in the center of the inlet 
flow channel was approximately 210 ft/sec. 

In the initial phase of the experiment tests were made at diffusor 
angles varying from /3 = 10° to/3 = 60°. At each angle, the flow rate 
was varied between the maximum of 2.4 and a minimum of 0.8 Ib^sec. 

Tests were made at each of the angles with smooth walls and with 
ribbed walls configured as indicated below: 

N: number of grooves = 10 

b: width of grooves = 0.125 in. 

t: spacing between grooves = 0.060 in. 

a: depth of grooves = 0.700 in. 

This rib configuration was chosen as a close approximation to the rib 
sizes used by Migay in the two dimensional diffusor study cited earlier 
and summarized in Table I. 

The results of this phase of the experiment indicated that the 
ribbed walls performed best, relative to the smooth walls, at ^ = 32 ^ 
and a velocity of approximately 155 ft/sec. Consequently, further 
tests were conducted at this operating point in an attempt to determine 
how the performance varied as a function of the number of grooves. In 
this phrse of the experiment the number of grooves was increased by 
one for each successive run. The results of these tests are discussed 
in Chapter IV. 

The program used to perform the required calculations is discussed 
in detail in Appendix A. A sample data sheet, a list of the variable 
names used on the sample data sheet, and a listing of input variables 
are included with the discussion of the program. 
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CHAPTER III 



COMPUTATIONS AND ANALYSIS 



3*1 Temperature measurement conversions ^ The thermocouple 
readings, recorded in millivolts, were converted to degrees Fahrenheit 
by the following equation: 



This equation is a least squares fit of the values for conversion 
from millivolts to degrees Fahrenheit as specified in the National 
Bureau of Standards Circular 561 for iron constantan thermocouples. 



at a temperature of 32 degrees Fahrenheit. 

3.2 Pressure measurement conversions . The barometric pressure, 
measured and recorded in inches of mercury, was converted to pounds per 
square inch by using the following equation: 



The two coefficients multiplying t, the temperature in degrees centi- 
grade, are thermal expansion coefficients for the mercury in the baro- 
meter and the brass upon which the scale of the barometer is engraved. 
The factor .4912 is the conversion factor for converting pressure 
measured in units of inches of mercury to units of lbs. per sq. in. 
Equation (3-1 the conversion factor, and the correction coefficients 
were obtained from Ref. 5, pages 1-157. 

3*3 F low rate measurements . References 6 and 7 provide general 
methods, procedures, and calculations concerning flow rate measurements 



32.t35-.98(rv,v')-.455'(wv')^ 



(3-1) 



Equation (3-1) is valid only if the reference junction is maintained 




(3-2) 
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using flow metering orifices and nozzles. The flow metering nozzle 
used to determine the flow rate for these tests was constructed in 
accordance with the specifications of Ref. 8. A nozzle was selected, 
rather than a square-edged orifice, because the calculations were felt 
to be more easily adaptable to computerized methods. 

From Ref. 8, page 5, Eq. (la), slightly modified, the flow rate 
through the flow metering nozzle in Ib^sec is: 

(3-3) 

The quantities in Eq. (3-3) are identified below: 

Vs/ flow rate (lbn,/sec) 

discharge coefficient 
^ density coefficient 

area of nozzle opening (in.^) 
a rea multiplier 

•(y Reynolds number correction 

conversion factor (32.17 ft. Ib^/sec^lbf) 
pressure drop across nozzle (psi) 
upstream "density" (non-standard units) 

By assuming the flow to be a perfect gas, ^ can be expressed as 
a function of the temperature, Tu. 

PvA 



6 .- 



RcTv 






(3-4) 



gas constant for air (53.35 ~ ^ ) 






temperature upstream (°R) 

The "equivalent flow rate," denoted as , and defined by Eq. 
(3-5), provided a computational means of determining the flow rate. 
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(3-5) 




In Eq. (3-5)>''<V has the units of in,^. 

After substituting for W and Eq. (3-5) becomes: 

W = « £ Ai-ft'i-'fT' 

' ' (3-6) 

The table below lists the values of the coefficients in Eq# (3-6) 
which are applicable to the particular nozzle installed for these tests. 



TABLE II 

SUMMARY OF FLOW RATE CALCULATION COEFFICIENTS 



Coefficient Symbol 


Coef f i cient Value 


Source of Value 


o< 


1,114 


Ref . 4 


€ 


Eq. (3-7) 


Ref. 2 


ft 


Eq. (3-8) 


Ref. 2 


■fv' 


1.0 


Ref. 4 




7.1728 


IT r’- 



P ^ Pa 

If X - f then the coefficient £1 can be expressed by the 



equat i on, 



1 



Z.^33 






i.Azs 



1-,3050-xV 






\AZ8 



(3-7) 



If is the temperature of the flow upstream of the 
pressed in degrees Fahrenheit, then the coefficient -f^ is 
equat i on, 

r , . _ _ /t«- foO\ 
ft . \+ .oo)b ( -i-oo ) 



nozzle, ex- 
given by the 



(3-8) 
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Equation (3-6) may be greatly simplified by combining all of the 
known constant coefficients into one constant, c. 



C - ®< -Pr Aci =• \ \ • 30 \r\}' 



Thus, Eq. (3-6) becomes: 

11.30 £-ft 



Pu. 



(3-9) 



Finally, the mass flow rate can be calculated by; 









(3-10) 



3.4 Flow conditions at the diffusor inlet . Velocity and total 
pressure measurements across the inlet of the diffusor indicated that 
the static pressure was uniform throughout the inlet plane. It was 
assumed that the total pressure at the inlet was constant in the region 
outside of the boundary layer. Based upon this assumption, the Mach 
number in the central region of the flow at the inlet was determined 
by solving the following isentropic relationship for ; 






The measured values of and were used to determine M 



(3-11) 



(3-12) 



The static temperature at the inlet was determined by using the 
previously determined Mach number and the Following equation: 



Ttv ^ 

T, 






T-1 






(3-13) 
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Solving Eq. (3-13) for yields the required relationship; 



T, - 






1 + 



(3-l‘>) 

Once Tj y the static temperature, was determined, it was possible 
to evaluate a, the speed of sound, and subsequently the midstream velo- 
city at the inlet. 






(3-15) 



The velocity in the central region of the flow is indicated by 

By assuming the flow to be a perfect gas, enough information was 
available to calculate the density of the flow at the inlet. It was 
assumed that the density remained constant across the inlet plane. 

p ; p. ^ 

' ( 3 - 16 ) 

Equation (3-16) expresses the density in units of slugs/ft^. 

The Reynolds number based upon the equivalent diameter of the 
inlet was determined by expressing the kinematic viscosity of air as 
a function of temperature. The expression represents a fit of the 
data tabulated on page 8 of Ref. 9» The following expression is con 
sidered valid only for the temperature range of 50°F. to 1 04®F : 

(o.fc t,+ UO)k 

(3-17) 

The units of are ft^/sec. 

3.5 Integrated f 1 ow rate at the di f f usor inlet . All of the cal- 
culations related to the conditions at the inlet and the performance 
of the diffusor were dependent upon the total pressure, static pressure, 
and total temperature measurements at the diffusor inlet. It was felt 
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that some method of evaluating the reliability of these measurements 
should be incorporated into the computational procedures. The method 
decided upon was to compare the flow rate at the inlet, computed from 
the velocity and density as determined from these measurements, against 
the flow rate measured at the flow nozzle. Since the two flow rates 
were determined from completely independent measurements, it was felt 
that any erroneous readings, or errors in recording the data, would 
become apparent when the two flow rates were compared. 

Prior to performing the calculations to determine the flow rate 
through the diffusor, tests were made with a hot wire anemometer to 
determine the velocity profiles at the inlet. Figure 6 is the velocity 
ratio, V/Vfnax^ plotted as a function of the distance, h, from the sur- 
face of the wall, Vpiax ^ velocity in the central region of the 
flow. Figure 7 represents the same velocity ratio as measured from the 
side wall. The profiles represented in Figs. 6 and 7 were measured at 
the mid-span position of the top and side walls respectively. 

In each instance, it was possible to find analytic expressions 
which closely approximated the velocity profiles. The data plotted 
were taken at several velocities representative of all anticipated 
operating conditions. Although it was recognized that the velocity 
profiles vary as a function of velocity, tests indicated that the vari- 
ation was minor and under the anticipated velocity operating range 
would not exceed the variations caused solely by instrumentation error. 
The total change in the flow rate determination over the extremes in 
the velocity profile predictions was found to be a fraction of 1 per 
cent of the total flow rate. 
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DIFFUSOR INLET VELOCITY PROFILE AT MID-SPAN OF TOP WALL WITH EXPONENTIAL 
EXPRESSION USED TO APPROXIMATE THE SHAPE OF PROFILE 
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DIFFUSOR INLET VELOCITY PROFILE AT MID-SPAN OF SIDE WALL WITH EXPONENTIAL 
EXPRESSION USED TO APPROXIMATE THE SHAPE OF PROFILE 



The expression used to approximate the velocity along the top wall 



is: 



rifcoci- yvMMV 



h = 1. 2.0 e 

( 3 - 18 ) 

When solved for V as a function of h, the distance from the wall, 

Eq. (3-18) becomes; 

l.Z0\^ 

(3-19) 

A similar expression was obtained for the shape of the velocity 
profile at the side wall* 



V = ^'-60 ] 



h = 0.60 e 






( 3 - 20 ) 



Solving for V as a function of h yields: 

Figure 8 is a cross-sectional view of the inlet plane. The region 
labeled "C" is the area outside of the effect of the boundary layer. 

It was assumed that the velocity in this region, previously referred to 
V_,^, was constant throughout. Measurements indicated that this 

m3X 

assumption was reasonable. 

A turbulence survey was made in the central region by feeding the 
hot wire anemometer output to an oscilloscope. The results of this 
study indicated that the 1 maximum velocity perturbation was less than 
1 per cent of the mean velocity, ^max* 

The regions labeled and “D” in Fig. 8 represent those regions 

within the boundary layer along the top, bottom, and side walls 
respectively. 
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INLET PLANE OF DIFFUSOR SHOWING THE LOCATION OF BOUNDARY LAYERS AND FLOW REGIONS 
AS USED IN THE DETERMINATION OF THE FLOW RATE THROUGH THE INLET 



Although it was not possible to accurately measure the velocity 
profiles on the bottom and opposite side wall, those measurements taken 
indicated that the variations from wall to wall and top to bottom were 
slight. Thus, it was assumed that the velocity profiles on opposite 
walls were identical. 

The corners, indicated as "E” in Fig. 8, were assumed to have a 
velocity distribution equal to the product of the ratios from 

each wall as a function of the respective distances. Figure 9 shows 
how this was done. The velocity assumed to exist at the point p, 
indicated U» Fig. 9, was assumed to be the product of the velocity 
ratio given by Eq. (3-19) for the distance hp from the top 

wall and the velocity ratio given by Eq. (3-20), for the 

specified distance, h 2 , from the side wall. 

The diagonal line indicated as “corner survey path line'' in 

# 

Fig. 9 is the line upon which a survey of velocities was made. Figure 
10 shows the comparison between what the above assumption predicted for 
the velocities along the line and the actual measured values of 
The analytic expression for the velocity at any point in the corner 
region is given as a function of the two distances h^ and h 2 « 

The purpose of this investigation into the velocity profile at the 
inlet was to obtain a measure of the reliability, or validity, of the 
pressure readings. This was accomplished by using the measured values 
of the velocity and density together with the estimates of the velocity 
profile to integrate across the inlet and determine the flow rate. 
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TOP OF DIFFUSOR 
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FIGURE 9 

THE CORNER REGION (E) OF THE INLET PLANE 
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The calculation for the mass flow rate is expressed by the follow- 
ing equation when the flow is considered to be incompressible: 

0 

Table III is a breakdown of the mass flow rate calculations for 
each of the regions specified in Fig. 9. 



TABLE III 

SUMMARY OF CONTINUITY CALCULATIONS 



Regi on 


Area 


Velocity 


Mass Flow (Ibj^sec) 


Eq. 


c 


0.008 


'^max 


? 3c A 


(3-24) 












D 


0.013 


Eq. (3-19) 


0.022 ^9c\ Y, =1^ 


(3-25) 












E 


0.020 


Eq. (3-22) 


0,028 f 5 


(3-26) 


F 


0.080 


Eq. (3-21) 


0-001 fScS'N*. 


(3-27) 



The flow rate determined by this integration was compared to the 
flow rate measured at the flow nozzle as discussed in Section 3»2. 

In general, the agreement was within 1 per cent, 

3*6 Diffusor performance . The conditions which existed at the 
diffusor inlet were found to be generally steady, and the measurements 
were felt to be reliable indicators of the flow at that plane. It was 
found, however, that conditions at the exit plane were, under most 
operating conditions, extremely turbulent. For large opening angles, 
of particular interest in this study, the flow actually reversed direc- 
tion along the diffusor walls after separation had occurred. Moreover, 
at opening angles greater than 26 degrees, the flow was unsteady in that 
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it would separate first on one side of the diffusor and then on the 
other. It was noted that under certain conditions the flow could be 



directed along one wall of the diffusor while complete separation ex- 



isted along the opposite wall. The flow could then be shifted to the 
opposite wall causing the separation to shift sides. It was also noted 
that the diffusor could be made to separate at both sides simultaneously 
causing a "free jet" flow straight out of the inlet throat. 

Because of these non-uniform conditions which existed at the exit, 
it was not possible to measure a single meaningful exit velocity. The 
only flow condition which could be accurately measured at the exit was 
the static pressure, Because the diffusor exhausted to atmospheric 

conditions at subsonic velocities, the exit static pressure was taken 
to be equal to atmospheric pressure. 

Not being able to measure an exit velocity posed a serious problem 
in attempting to measure diffusor performance. Most definitions of 
diffusor efficiency, for example, require a knowledge of the exit velo- 
city, \l2m Migay, Ref, 1, defines diffusor efficiency as: 



He makes no attempt to explain the meaning of V 2 or how it was obtained, 
A more convenient and perhaps realistic criteria for diffusor per- 
formance are the coefficient of pressure recovery, defined as 



This measure of diffusor performance has the distinct advantage 
that all of the quantities can be easily determined, velo- 

city of the flow in "Region C," Fig, 8, was chosen as the reference 




(3-28) 



C 



2 ^ 





(3-29) 



velocity V], to be used in Eq. (3-29) 
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In order to compare the results of this test with those published 
in Ref. 1, some scheme had to be decided upon to produce a quantity V 2 
for Eq. (3-28). One method used was simply to assume that the central 
flow at the inlet, with velocity was expanded completely and in- 

compressibly from area A] to area A 2 * ^2 ideal expansion of the 

midstream velocity at the inlet. 

A second means of evaluating the efficiency was to determine the 
ratio of the measured to the ideal coefficient of pressure recovery. 

The ideal coefficient of pressure recovery can be expressed as a 
function of the inlet and exit areas, 






(3-30) 

Equation (3-30) was derived by substituting Bernouilli's equation 
and the equation of continuity for incompressible flow into Eq. (3-29). 

3.7 Measurement re 1 i abi 1 i ty i ndi cator . In order to obtain an- 
other indication of the reliability of the pressure and temperature 
measurements at the diffusor inlet, a comparison was made of the mea- 
sured flow rate based upon the flow metering nozzle and the flow rate 
calculated by integrating the velocity over the area of the diffusor 
inlet as described in Section 3.^. The quantity referred to as the 
’’data reliability indicator" was defined by the equation: 



w. 



:: \ - 



(3-31) 
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CHAPTER IV 



RESULTS AND RECOMMENDATIONS 

Results . This project was conducted in three phases. The 
pre 1 imi nary phase consisted of planning, building and evaluating the 
test equipment. The second phase consisted of the development of 
procedures and data reduction calculations. The final phase encom- 
passed the testing required to determine the most desirable operating 
range for conducting further tests and an investigation to determine 
the effect of the number of grooves upon the ribbed wall diffusor per- 
formance. The results of each of these phases of the project are dis- 
cussed individually in the following sections. 

4.1.1 Preliminary phase . The comparison of the flow rate meas- 
ured at the nozzle and the flow rate computed at the diffusor inlet 
provided a cross check of the pressure readings and the computed 
values of the Mach number and velocity. The data reliability indica- 
tor defined by Eq. (3-31) reflected the comparison of these two flow 
rates. The average value of this indicator for the last 50 data 
points recorded during the project v\fas 0.0044. This indicates that 
the agreement was, on the average, within 0.5 per cent. This is com- 
parable with the accuracy of the flow nozzle which is - 0.5 per cent. 

Section 3*6 of chapter III discusses the two methods used to 
determine the diffusor efficiency. It can be shown that the ratio 

is identical to the efficiency defined by Eq. (3-28). Both 
methods of determining the efficiency were included in the calculations 
dealing with diffusor performance since Eq. (3-28) is based upon the 
calculated value of the density,^ , whereas the ratio of pressure 
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recovery coefficients does not depend upon . The comparison of 
these two efficiencies provided, therefore, an indication of the 
reliability of the measured value of the temperature and pressures at 
the inlet. 

The third indication of the reliability of the test equipment 
and procedures was the repeatability of the tests. Table IV shows 
the results of tests run on separate days under similar operating 
conditions. It was felt that the repeatability indicated in Table IV 
was reasonably good and that a realistic performance determination 
could be made on the basis of these tests. 

The results of the tests seemed to indicate that the test equip- 
ment and procedures could be relied upon to an extent consistent 
with the intended use of the calculated results. It was felt that 
sufficient safeguards had been incorporated into the computational 
system to provide adequate warning indications in the event that 
erroneous data were being processed. 

4.1.2 Determination of the most desirable operating range . To 
accomplish the objective of this phase of the study, ribbed wall and 
smooth wall tests were run over a wide range of operating conditions. 

Table V summarizes the results of this phase of the tests. Al- 
though runs were made from = 10^ to ^ = 60 ° for velocities 
ranging from 65 ft/sec to 210 ft/sec, only those values of ^ which 
were of particular interest to this study are presented in Table IV. 

Figure 11 illustrates the operating range of particular interest 
for ribbed wall investigations. Figure 12 shows how the ribbed wall 
efficiency compares to that of the smooth wall diffusor in this region. 
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TABLE IV 



RESULTS OF REPETITIVE TESTS UNDER SIMILAR 
OPERATING CONDITIONS 



n 


Ribbed Wa 1 1 


Smooth Wa 1 1 


Run No. 




Run No. 


7s 


26° 


219 


61.75 


220 


65.30 


26° 


221 


62.86 


202 


66.53 


O 

00 

C\J 


219 


59.63 


220 


58.58 


o 

00 

CNJ 


221 


59.68 


— 




o 

o 


219 


58.94 


220 


53.51 


o 

o 


221 


59.17 


204 


53.60 


32° 


219 


56.96 


220 


47.84 


32° 


221 


57.36 


207 


47.17 
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table V 



COMPARISON OF SMOOTH WALL AND RIBBED WALL 
DIFFUSOR EFFICIENCIES* 



/3 


7s 






20 ° 


83.06 






22 ° 


81.21 


78.29 


0.9640 


24° 




69.79 




26° 


65.30 


61.75 


0.9456 


28° 


58.58 


59.63 


1.0179 


O 

o 


53.51 


58.44 


1.1015 


32° 


47.84 


56.95 


1.1906 


34° 


43.51 


34.55 


0.7941 


36° 


39.30 


29.82 


0.7588 


38° 


34.55 






40° 


31.66 


27.69 


0.8746 



*A11 data recorded at a measured flow rate of 
approximately 1.78 Ib^/sec. Ribbed wall con- 
f i gurati on; 

A = 0,7 in. 

B = 0.125 in. 

TH = 0.06 
N = 10 



42 




1 1 1 1 1 1 ► 



a^ 


C30 




vO 


• 


• 


• 


• 


o 


o 


o 


o 



u^ 

• • • 

o o o 






43 



SUMMARY OF RIBBED WALL AND SMOOTH WALL DIFFUSOR EFFICIENCIES 



table \l 



COMPARISON OF SMOOTH WALL AND RIBBED WALL 
DIFFUSOR EFFICIENCIES* 



/3 


Vs 


Vr 




20° 


83.06 






22° 


81.21 


78.29 


0.9640 


24° 




69.79 




26° 


65.30 


61.75 


0.9456 


28° 


58.58 


59.63 


1.0179 


O 

o 


53.51 


58.44 


1.1015 


32° 


47.84 


56.95 


1 . 1906 


34° 


43.51 


34.55 


0.7941 


36° 


39.30 


29.82 


0.7588 


38° 


34.55 






40° 


31.66 


27.69 


0.8746 



*A11 data recorded at a measured flow rate of 
approximately 1.78 Ib^/sec. Ribbed wall con- 
f i gurati on: 

A = 0,7 in. 

B = 0.125 in. 

TH = 0.06 
N = 10 
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SUMMARY OF RIBBED WALL AND SMOOTH WALL DIFFUSOR EFFICIENCIES 




CN O On 00 

• • • • • 

^ o o 
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RATIO AS A FUNCTION OF /3 TO COMPARE THE RIBBED WALL DIFFUSOR 

PERFORMANCE TO THAT OF A SMOOTH WALL DIFFUSOR 



4.1.3 Effect of vary 1 ng the number of grooves upon r i bbed wa 1 1 
di f fusor performance . Table VI is the result of the third and final 
phase of this investigation. Unfortunately, time did not permit this 
phase of the testing to proceed beyond N = 10 grooves. One interesting 



fusor decreased by approximately 20 per cent as the number of grooves 



Migay, Ref. 1, reports a similar drop in efficiency for a three di- 
mensional conical diffusor operating at an opening angle of 31 degrees. 
His explanation for this decrease in efficiency is that it is caused 
by "acoustic irradiation of inner-fin hollows." No attempt was made in 
this study to either explain or further investigate the phenomenon. 

4.2 Recommendat i ons . Much experimental work remains to be done 
in the investigation of ribbed walled diffusors. It is hoped that 
continued studies in this field will make available data upon which a 
theoretical analysis of ribbed wall diffusor performance can be based. 

It is felt that the next step in the investigation of ribbed 
wall diffusors should be an attempt to obtain data reflecting the 
effect on performance as other parameters of the rib configuration 
are varied. 

In conjunction with further experimental studies into the effect 
of ribs on the performance of diffusors, it is felt that studies into 
the mechanism of flow past grooves or ribs should be thoroughly in- 
vestigated. Fox, Ref. 10, has done much work in this field. It is 
felt that the results of his study would be of great value in attempt- 
ing to understand the mechanism of the ribbed wall diffusor operation. 



result of these tests is illustrated in Table VI. 





noted that the efficiency of the ribbed wall dif- 



was increased from 7 to 8. This phenomenon did not occur 
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TABLE VI 



EFFECT OF VARYING THE NUMBER OF GROOVES UPON THE PERFORMANCE 
OF A RIBBED WALL DIFFUSOR.* 



fi - 28°, A = 0.7, B = 0.125, TH = 0.06 



%/ y )^ 


0.7979 


0.9339 


0.9507 


0.9706 


1.0179 


N 


6 


7 


8 


9 


10 



- 30°, A = 0.7, B - 0.125, TH = 0.06 



?/?/% 


0.7443 


0.9587 


0 *8866 


1.0525 


1.1015 


N 


6 


7 


8 


9 


10 



= 32°, A = 0.7, B = 0.125, TH = 0.06 





0.5847 


0.9505 


0.7306 


1.1378 


1.1906 


N 


6 


7 


8 


9 


10 



*All runs for this test made with measured flow rate of 
approximately 1.78 Ib^/sec. 
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The possible use of ribs, or grooves, as a means of affecting 
boundary layer control in a variety of external as well as internal 
flow situations seems to indicate that further research into the 
behavior of ribs is warranted. It is hoped that the equipment, pro- 
cedures, and experience assembled in the course of this study will 
serve some useful purpose in any such subsequent investigations. 
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APPENDIX A 



COMPUTER PROGRAM 

Genera 1 statements concerni ng program " Di f fusor All of the 
calculations discussed in Chapter III were programmed for the IBM 
36 O -67 Digital Computer in operation at the Computer Facility, Naval 
Postgraduate School, Monterey, The program was written in Fortran IV 
language* All of the calculations as discussed in Chapter III are 
utilized in the program named "Diffusor*" Appendix B is a listing of 
the program. Included in the program listing are "comment" cards to 
help explain each step of the calculation. A basic understanding of 
the Fortran IV language should be a 1 1 that is required to follow the 
program operation. In this Appendix the emphasis is directed toward 
the input and output phase of the program operation. 

Program input . The special listing of symbols entitled "Program 
Input Variables" included with this Appendix defines the Fortran 
variables names used as input for program "Diffusor." Figure 13 is a 
sample data page illustrating the order and variable names as recorded 
during the tests. Figure 14 shows the field specification for each of 
the input variables. 

The first data card of each program run is a 72 space "Hollerith" 
statement. The first 72 spaces of this card may be utilized to include 
supplementary or amplifying statements to be printed as part of the 
program output. Any information punched on this card will be repro- 
duced directly beneath the title line on the output page. 



49 



PROGRAM INPUT VARIABLES 



MM 


Number of sets of data to be read 


IRUN 


Number of run 


MO 


Month (numeric) 


I DAY 


Day of month (numeric) 


PATM 


Atmospher i c pressure 


PINOZ 


Static pressure ahead of nozzle ( ) 


OLTP 


Pressure drop across nozzle (AP ) 


TN 


Total temperature at nozzle 


TEMP 


Temperature (*^C) for PATM correction 


PI 


Static pressure at inlet P^ 


ai 


Dynamic head at inlet (P^-Pj) 


TTI 


Total temperature at inlet ( 


BETA 


Total diffusor opening angle 


TPR 


Outside ambient air temperature (‘^F) 


AA 


Depth of grooves 


BB 


Width of grooves 


TH 


Spacing between grooves 


NN 


Number of grooves 
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TURBO-MACHINE data PGS 145 (4-66 ) 



I 1J_ 

SAMPLE INPUT/OUTPUT DATA FOR PROGRAM DIFFUSOR 



72 Space Hollerith Statement Card 



1 2 




1 




MM 





1 3 


4 6 


7 9 


10 19 


20 29 


30 39 


40 49 


50 59 


207 


02 


13 


29.83 


12.35 


18.95 


1.50 


22.0 


IRUM 


MO 


IDAY 


PATM 


PINOZ 


DLTP 


TN 


TEMP 



1 10 


11 20 


0 

CM 


31 40 


41 50 


- 0.850 


2.155 


1 .50 


32.0 


57.5 


PI 


Q1 


TTl 


BETA 


TPR 



1 10 


1 1 20 


21 30 


31 40 


0.70 


0.125 


0.06 


10 


AA 


BB 


TH 


NN 



FIGURE 14 

INPUT DATA FIELD AND FORMAT SPECIFICATIONS FOR PROGRAM DIFFUSOR 
(Top line indicates card spacing) 
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The second data card of each computer run will indicate the 
number of data sets which are to be read* The name given to this 
"fixed point" variable is MM* 

The remaining data cards represent the input to the program* 

Each data point requires three cards* As illustrated in Fig* l4, the 
first of each series of three input cards contains the following: 

IRUN, MO, IDAY, PATM, PINOZ, DLTP, TN, and TEMP. The second of each 
series of three cards includes: PI, QI, TTI, BETA, and TPR, The last 

in the set of three cards used for each data point contains the fol- 
lowing information related to the ribs of the diffusor: AA, BB, TH, 

and NN* The symbols used to identify the input variables are defined 
in this Appendix* For smooth wall runs, this card was left blank 
resulting in zeros for all of the rib parameters* 

Special attention is directed to the fact that the "I" format 
variables are right justified within their specified fields whereas 
the "F format variables can be punched anywhere within the specified 
field provided the decimal is punched with the number* The field 
specifications and the format for each of the input variables is 
included in Fig* 14* 

Program output * Figure 15 is a sample of the output from program 
"Diffusor*" The same data used to illustrate the input in Figs* 13 
and l4 are printed as output in Fig* 15. 

The second line of output is the Hollerith statement read into 
the program as the first data card* 

The grouping labeled DATA CHECK is a print out of the raw data 
as read into the program* It was used to check the output for 
possible errors in card punching* 



53 



sample input/output data for program niFFUSDR 



u 





GC 








LAC LO 












V 




1 








H w cr 




r' 








U LA -J 




•-H 












1 

tv 








c c 










4 


clc c occfn>tr^r- 








o *o 


0 


ir 1 r> ^ 




LL 




azLL 


r- 


fvLLU c u^c cc <4 






u 




• 


•oc oc •••••• 




<J 




V 


c 


a >4 ore cN'^c c c 




c 


c 


OCOOL 




oc f\,C 






u 












c 






• • 










II 


II II II II II II II II II 








C C rr 










c 


^ ^4 




• ^ 






c 
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^<1^ •CM 






• 


c cc 
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OU_> 








• • • 


LL* 


^ QC c OC H- te ^ 






CsJ 


•— i(NJ^ 


> 


OC •D— <IC-J-J 


►- 




II 




c 


a o>OKa u<i<i 


<I 




a 


II II II 


o 


O O LU lULU 


c_ 




X 


u. 


LU 


O JLUK-J3>CC 






u. 


K 


CX 


<i^ujzca^~^ 


z 




h 


<I 




Jar j>— jUj 


c 






or LU 


uj**<k •2^u->oo: 








^ ►- 


oTL/^azc^K oza 


►— 








^ z uj z ^ OJ ^ 








LoCa 


t/^CZO >tt/'U-i/>QC 








< j 
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LUwtr •LL<ia:Z^a 


Zi 
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XU.3 


U,KO>CO^X ww 


-a 




IT 


c 


CXi— 


W.KJZ • 


> 
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K J 
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— O 


LL 






UJ2LL 


ZX i>nz X»-'U-ii.-eo 








KU. 


u-C<Z > o»->-a:zz 


UJ 
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<1-JC 


COXU-0<t<I<ia LL LL 


O 




z 


QC <IU. 




>CaXa:«j 


z 




I- 


>0L 


•h“ 


joc 


< 








LLUJ 




3- 






oru- 


U--J 


1 LL^OOLL 


O 




c 


-JCO 


CZ 


Qc-aCiLU 


C 






• •LLLL'QC 




<0»-'LU0 



u 

oc 



a 



o 



3 

a 



o • • • 

C ^ 



CJC c c 

c <. 

• • • • 

CNJOC c c 



II II II II II 



Of^ 

irnr 

a II 

•a 

oca 



II 

a 



JO 

Cir. 



fSN 

^si 

c 



3 

c 



a ac 



<r<i< 

CO^i— *— 

a v/Jt/;cr 

o a a 



<i 

or 

c 



c 



r 

i- 

c 

< 

e 

II 



< 

O' 



< 

or 

3 

C 



01 

C3 







CK 


J 
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— «• 




ITK 
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c irT-c 
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l~ 
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SAMPLE OUTPUT PROGRAM DIFFUSOR 



The input sections under the headings of "Flow Rate" and 
"Diffusor Performance" are the converted values of the input data. 

The output sections under these two headings are those values 
calculated in accordance with the discussion in Chapter III. 
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APPENDIX B 



FORTRAN IV f. LEVFL 1. MOO 0 MAIN OATF = 68061 



C 

C ASSIGN VALUES TO FIXEO CONSTANTS NEEOFO FOR COMPUTATIONS 
C 



OOCl 




TST0=5ie.7 






PST0=2^.92^. A921 


00^3 




TEXAS =.8 


COCA 




32.17 


CC0*5 




53.35 


V C C 6 




GAM=r 1.4 


CCC7 




MM = 0 


OOCfl 


c 


A=2.*G*778.3 




C 

r 


READ IN INPUT DATA AND PRINT OUT FOR DATA CHECK 


0CC9 


L 


READ CStlO^O) 


CC1<^ 




READ IS.l'') NUM 


0011 




15 REAP C5f2C) TRUN, MO, I DAY, PATH. P1N07 . DLTP . TN. TEMP 


0^12 




READ (5,251 P 1 , Q 1 , TT 1 , BETA , TPR 


0 013 




READ (*^,26) AA,RR,TH,NN 


OOIA 




MM=^MM4>1 


C015 




WRITE (<S,3<3Q) 


c:i6 




WRITE (6,A'?0) 


CC17 




WRITE (6,1^001 


C Cl 8 




WRITE C6.A01) IRUN,MO,IDAY 






WRITE (6,A^2) 


0020 




WRITE C6,A03I BETA 


0‘*^21 




WRITE C6,A''A) PB 


0022 




WRITE (6,A'^S)AA 


0023 




WRITE (6,A(‘6) TH 


002A 




WRITE C6fAr.7) NN 


G'>25 




WRITE C6,3B0) 


r026 




WRITE CB,3R1) P ATM , P 1 NOZ , OLTP , TN , TEM P 


J 0 27 


r 


WRITE (6,382) PI ,01 ,TTl,TPR 




C 

r 


CONVERT DATA INTO DIMENSIONALLY CORRECT UNITS HE MEASURE 


00 2 8 




PATM=PATM>^( l,-( (,0Cni634*TEMP) /( l. + .0^'^1818*TEMP ) ))*, A^l 2 


0029 




o2=PATM 


003^ 




Dl TP=0LTP*TEXAS*.03613 


00 31 




R1N0Z = PATM4-P1N0Z*TEXAS*.03613 


C022 




01=Q1*TEXAS*. 03613 


0033 




P1-PATm^pi*.03M3*tFXAS 


C03A 




PT1 = P1'^01 


003*^ 




TT 1= 32. ♦35.98*TT1-. A35*TTI**2^A59, 7 


C036 




TN= 32.^3S.9R*TN-. A35*TN*^2 


CC37 


r 


TNR = TN«-A59. 7 




c 

c 


BEGIN FLOW RATE CALCULATIONS USING FLOW NOZZLE READINGS 


c :3P 


X=nLTP/PlN07 


C039 




Y=( 1.-X)**i.42P 


GOAu 




Z=( l.-X)**0.?86 


OCAl 




EPS=S0RT(2.A33*Y/X*( 1.-7) /( 1.-. 305*Y) ) 


00A2 




E T= !.♦. Geo'll 5* (TN-6C. ) 


00A3 




WEOI V=ll. 3*EPS*ET«S0RT( X) 


OCAA 


r 


W=WE0IV*PlN074cS0RT(G/(R4tTNR) ) 




C 

r 


COMPUTE MACH NUMPER 


00A5 


v> 

r 


AMACH=S0RT(2./(GAM-1. ) ♦( ( PTl/Pl )♦♦( (GAM-1. ) /GAMI-l. ) ) 




c 

r 


CCMPUTE STATIC TEMPERATURE AT INLET 


00A6 


r 


T1 = TT1/ ( 1. -*.(GAM-l. ) /2.*AMACH**2) 




c 

r 


CCMPUTE VELOCITY IN THE CENTRAL FLOW REGION 


C0A7 


L- 

r 


VMAX = AMACH*SORT(GAM»R*G>*tTl ) 




C 

r 


COMPUTE THE DENSITY OF THE FLOW IN THE CENTRAL REGION 


0CA8 




R0 = PI*144. /(R*G*T1 ) 



C 



05 / 12 / 3 '' 
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FTRTRAN IV G LEVEL U MOO 'I MAIN 



OATF ^ ''5/l?/?0 



C050 

0051 

C052 

0C53 

0054 

0055 

0056 

0057 
0^58 
0C59 
C06? 
0C61 
C062 

0063 

0064 

0065 

c :66 

0'^67 

0C68 

006 ^ 

^ 07 (. 

t^C 71 



0072 



OJ73 
03 74 
0075 
C376 



0277 

0378 



0 C 75 



0C82 

CC81 



0 :'i 82 



0383 



0 184 



0 J8 5 



C BEGIN THE NtIMFRICAL INTEGRATION TO DETERMINE THE FLOW PATE AT TWP 
C iNLET PLANE 
C 

W$UBD=O.D 

WSUBE= 

WSUBF= 0,0 
H =-.025 

00 ICO N=ltll 6 f 5 
h=H +.05 

V = VMAV*( l.-.ni 667 *C ALOGI 1 . 2 /Hn** 2 ) 
no W$UBF=Rn*G’(cV*,r 5 * 4 , 8 / 144 . 4 WSURF 

WF=2. ♦WSUBF 

WC=R 0 *G*VMAX^ 4 , B+ 1 . 6 / 144 . 

H =-.025 

on 2 C 0 N=l, 56*5 
H=H +,05 

V=VMAX*C 1 ,-. 03333 *( ALOGI. 60 /H )!♦♦?) 

WSUBD=RO*G«V*.r 5 * 1 . 6 / 144 . ♦WSURO 
Hl =-.025 

00 20 '' J = l. 116,5 
Hl=Hl +.05 

V = VMAX*( l.-.O 1667 *( ALOGC 1 . 2 /Hl) I **2 ) ♦ (1 . - n 33 3 3 ’*' f ALOGI , 6 /Hn ^*2 ) 
2)0 W$URE = Rn*G=^V*, 0025 / 144 .^ WSUBE 

WE=W$UBE« 4 . 

wn=?,*wsuRn 

WCALC=WC+Wn+WF^WF 

c 

C COMPUTE THE COEFFICIENT OF PRESSURE RECOVERY 
C 

CPR=:<P 2 -Pn /01 

C 

C CCMPUTF THE AREA RATIO OF THE DIFFfJSOR FOR THE SPECTFIEO OPENING 
C ANGLE 
C 

8 ETA=BETA/ 57.3 

ALFA=RETA/ 2 , 

RPTA=RETA* 57, 3 

APATI 0 =. 33333 *( 3 . 1 fi 75 + 16 . TANI ALFA)-. 1875 /COSCAL fa n 
C 

C COMPUTE REYNOLDS NUMBER BASED UPON MIDSTREAM VELOCITY AND THE 
C HYDRAULIC DIAMETER OF THE INLET FLOW CHANNEL 
C 

VI S = *n 00 ''' 0 ) 6 *( TT 1 - 45 Q. 7 ) + #00012 
RE=VMAX«, 4 /VIS 
C 

C COMPUTE the ideal COEFFICTENT of PRESStJRF RFCOVERV 
C 

PRSTAR= 1 .-( 1 ./ARATI 0)**2 

C 

C COMPUTE THE referred FLOW RATE (NGRMALI 7 ED TO STANDARD DAY 
C CONDITIONS) 

c 

TPR=TPR+45R, 7 

W$TD=W*SORT( TPP/TSTD) /(PATM/PSTO) 

C 

C COMPUTE the VELOCITY AT THE EXIT BASED UPON AN IDEAL EXPANSION 
C OF THE mid stream FLOW 
C 

v 2 =vmax/arat in 

c 

c CCMPUTF EFFICIENCY BASED UPON THE DEFINITION OF E 0 ( 3 - 28 ) 

C 

EENCY=( P2-oi)#144. /( .5*R0*C VMAX*#?-V2**2) ) 

C 

C COMPUTE EFFICIENCY BASED UPON THE lOFAL COEFFICIENT OF PRFSSUPF 
C RECOVERY 
C 

EFFN 2 = CPR/PRSTAR 

c 

C determine THE DATA RELIABILITY INDICATOR 
C 

ERR= l.n - CW/WCALC ) 

C 
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FCRTRAM IV G LFVFL It MOO 0 



MAIN 



OATF = 68061 



05/12/3^ 



CC86 

0087 

C088 

0C8S 

OC<50 

UC91 

00^2 

0093 

0094 
cn<;5 
0096 
C097 
0098 
0C99 
OIGO 
ClCl 
Olf 2 
0103 
01C4 
01C5 
0106 
0107 
01C8 
C1C9 
0110 
Gill 
0112 

0113 

0114 
0116 

0116 
0117 
Cl 18 
U119 
012C 
0121 
C) 122 

0123 

0124 
G12 6 

0126 

0127 

0128 
C129 

0130 

0131 

0132 

0133 

0134 

0135 

0136 

0137 

0138 

0139 
C140 

0141 

0142 

0143 



C 

C 



PRINT ALL VALUES AS COMPUTED 

WRITE (6»438) 

WRITE (6f409l 
WRITE I6»41DI PlNOZfW 
WRITE (6f411l DLTPfWEOIV 
WRITF (6f412) TNtWSTO 
WRITE (6,413l 
WRITF (6.409) 

WRITE (6.414) PATM.CPR 
WRITF (6.419) PI 
WRITF (6.416) PTl.VMAK 
WRITE (6.417) P2.RC 
WRITF (6.418) RF 
WRITF (6.419) AMACH 
WRITF (6.420) ARATIO 
WRITF (6.421 ) WCALC 
WRITE (6.423) PRSTAR 
WRITF (6.424) FFNCY 
WRITF (6.429) EPFN2 
WRITE (6.426) 
write (6.427) ERR 
IF (MM.LT.NUM) go to 15 

10 FORMAT ( I?) 

20 FORMAT (3I3.9F13.4) 

25 FORMAT (9F10.4) 

26 FORMAT (3F1C.4.I10) 

380 format (15X.'DATA CHECK'/) 

381 FORMAT (25X. •PATM=» .F6.3. • P 1 NO/ = • , F6 , 3 . • 

1 TN='.F5.2.' TEMP=' .F5. 2. • DFG C • ) 

382 FORMAT (25X,'P1 = %F6.3.' 01='.F6.3.* TT1 = «,F8,2.' 

15.2// ) 

399 FORMAT (IHl/////) 

400 format (30X^ 'TWO DIMENSIONAL RIRRED WALL niFFMSOR PERFORMANCE EV 
lALUATION DATA'//) 

401 format (15Xt4HPUN . I 3 . 80 X . 5HDAT E . I 2 . IH- . I 2 . 3H-68/ / / / / ) 

40? FORMAT ( 15X. 'CONFIGURATION: '/) 

403 FORMAT (25X.'BFTA (DIFFUSOR ANGLE) ='.Fir.4.' 

404 FORMAT (29X^'WIDTH OF GROOVES (B) ='.F10.4.' 

405 FORMAT (25X, 'DEPTH OF GROOVES (A) ='.F1^.4.' 

406 FORMAT (29X. 'SPACING BETWEEN GROOVFS (T) ='.Fir-,4.' 

407 FORMAT (?9X,»NUMBER OF GROOVES (N) = '.13//) 

408 FORMAT (15X.'FL0W RATE. MEASURED WITH 3.022 INCH FLOW NOZZLE'/) 

409 FORMAT ( 2C X I NPUT :'. 3^X .' OUTPUT :' ) 

410 format (25X.'P1(N07ZLE) ='.F8.4. 

ISS) =',F7.4.' LP/SEC) 

411 FORMAT (25X. 'DELTA PRESSURE ='.F8,4. 

ICW RATE=' .F7.4t ' SO. IN. ' ) 

412 format (29X, 'TOTAL TEMPERATURE ='.F8.4. 

1 RATE =',F7.4,' LR/SEC'//) 

413 FORMAT ( 1 5 X . ' 01 FFUSOR PERFORMANCE DATA'/) 

414 FORMAT(25X. 'ATMOSPHERIC PRES. =',F8.4. 

ISSURE RECOVERY =',F9.4) 

415 FORMAT (25X.'P1 (INLET STATIC) =^'.F8.4. 

HONS: ' ) 

416 FORMAT (25X,'PT1 (INLET TOTAL) ='.F8.4. 
lENTRAL CORE) =',F9.4.' FT/SEC) 

417 format (25X,'P? (EXIT STATIC) ='.F8.4. 

KCENTRAL CORE) ='.E9.4.' SLUGS/CU FT' 

418 FORMAT (76X. 'REYNOLDS NO. ( EOV. D T A. ) = ' . E9. 4 ) 

419 FORMAT (76X.'MACH NO. (INLET CORE) ='.F9.4) 

420 FORMAT (71X.'AFEA RATIO (EXIT/INLET) =i,f9.4) 

421 FORMAT (71 X. 'CALCULATED MASS FLOW RATE ='.F9.4.' 

423 FORMAT (71X. 'IDEAL PRES. RECOVERY COEF, =',f9.4) 

424 FORMAT ( 7 1 X . ' E FF I C 1 ENC Y (USING IDEAL V2) ='.F9.4) 

425 FORMAT ( 7 1 X . ' EFF I C I ENC Y (PR/PR IDEAL) ='.F9.4) 

426 FORMAT (//15X.'DATA RELIABILITY INDICATOR'/) 

427 FORMAT (25X.'1.0 - (FLOW RATE MEASURED/FLOW RATE CALCULATED) ='.F7 
1.4) 

1J30 FORMAT (72H 

1 ) 

END 



OLTP=' .F7,3.' 

TPR=' ^F 



DEG' ) 
IN. ' ) 
IN. ' ) 

IN. ' ) 



PSIA'. 13X. 'FLOW RATE ( MA 
PSIA '. 13X, ' equivalent FL 
DFG F' . 12X. 'REFER FD FI OW 

PSI '.lix. ' COEF. OF PRF 
PSIA'.12X.' INLET CONPIT 
PSIA'.12X,' VMAX (C 
PSI A' . 12X. ' PENS TTY 



LB/SFC ) 
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